A dynamical dark energy model with a given luminosity distance 
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It is assumed that the current cosmic acceleration is driven by a scalar field, the Lagrangian of 
which is a function of the kinetic term only, and that the luminosity distance is a given function of 
the red-shift. Upon comparison with Baryon Acoustic Oscillations (BAOs) and Cosmic Microwave 
Background (CMB) data the parameters of the models are determined, and then the time evolution 
of the scalar field is determined by the dynamics using the cosmological equations. We find that the 
solution is very different than the corresponding solution when the non-relativistic matter is ignored, 
and that the universe enters the acceleration era at larger red-shift compared to the standard ACDM 
model. 
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According to our current understanding about the cos- 
mos, we live in a flat universe which expands in an accel- 
erating rate and it is dominated by a dark component. 
Identifying the origin and nature of dark energy is one 
of the biggest challenges for modern cosmology. Since 
we can only speculate about what dark energy could be, 
many cosmological models have been proposed and stud- 
ied so far. The simplest candidate is a cosmological con- 
stant with state parameter w = —1, but models with 
a dynamical component with an evolving state param- 
eter w{z) also exist in the literature (for a review on 
dark energy models see e.g. [l|). Scalar field models with 
a non-canonical kinetic term have been discussed in k- 
inflation [4I, in which inflation is not due to the poten- 
tial, but rather to the kinetic term in the lagrangian, and 
in k-essence models. [3,|j|; which are designed to address 
the issue why the cosmic acceleration has recently begun. 
Very recently, in [5] the authors studied a k-essence model 
assuming i) a closed form parameterization for the lumi- 
nosity distance (for distances in cosmology see e.g. (y]), 
proposed previously in [7| , and ii) that the Lagrangian of 
the scalar field does not depend on the field itself. How- 
ever, in that work the matter has been ignored, and in 
the present article we wish to show that if the presence 
of matter is also taken into account, the time evolution 
of the scalar field is very different. We also find that 
according to this model, the universe enters the acceler- 
ation era at a larger red-shift compared to the standard 
ACDM model. 

The theoretical framework is basically the same as 
in [5|, apart from the inclusion of matter in the cosmo- 
logical equations, and the use of BAOs and CMB data 
to be mentioned later on. However, we start by giving 
a summary of all the ingredients we shall be using, so 
that the present work is self-contained. We focus here on 
the case of a fiat universe, since this is a robust predic- 
tion of infiation [^, ^ , and it is also supported by current 



data [lO[. The framework is based on four-dimensional 
General Relativity coupled to a single scalar field with a 
general Lagrangian £(0, X), where (p is the scalar field 
and X — (dcf))'^ /2 is the standard kinetic term. There- 
fore, our model is described by the action 



S = 



IBttG 



£(</), X) 



(1) 



The energy-momentum tensor for the scalar field is given 
by 



Tj^.t) = C.xd. 



^9tJ.u 



(2) 



where ,X denotes differentiation with respect to X. We 
can recast this energy-momentum tensor into the form of 
the energy-momentum tensor for a perfect fluid 



T^l-^'> = ip + pWu, 



P9ti:^ 



(3) 



where p, p are the energy density and the pressure of the 
fluid respectively. The hydrodynamical quantities p, p, u^ 
are given in terms of 0, X as follows 



p = C 

p = 2Xp^x - P 
dn(t> 



(4) 
(5) 

(6) 



Therefore, the equations of motion for the system grav- 
ity +scalar field arc just the first Friedmann equation and 
the equation for energy conservation 

SttG 



H' 



-iH{p + p) 



(7) 
(8) 



where H = a/a is the Hubble parameter, a is the scale 
factor and the overdot denotes differentiation with re- 
spect to cosmic time. Defining the sound speed of the 
scalar field 
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P.X 
P,X 



1 + 2A: 



P,xx 
P,x 



(9) 



the equation for energy conservation takes the form 
4> + 3clH<j> + -^ = 



P,x 



(10) 



which generahzes the usual Klein-Gordon equation of a 
canonical scalar field in a Friedmann-Robertson- Walker 
background. 

If we now include also the non-relativistic matter, the 
equations of motion for our system are the Friedmann 
equations and the conservation equation 



H' 



■ P 



H = --{p+p) 
= p + SH{p + p) 



(11) 

(12) 

(13) 

where k^ = SttG, G is Newtons constant, /?„ is the energy 
density of matter (pm = 0), px,Px are the energy density 
and pressure of dark energy respectively, p = pm + Px is 
the total energy density and p = Pm + Px = Px is the 
total pressure. Assuming no interaction between dark 
energy and non-relativistic matter, we have two seperate 
conservation equations 



= p„r + 'SHpjn 
0^px+iH{px+Px) 



(14) 
(15) 



where the second equation for dark energy is equivalent 
to the equation of motion of the scalar field. If we further 
assume that the Lagrangian of the scalar field depends 
on the kinetic term only, C = —F(X), the equation of 
motion can be integrated once 



XF^ 



Ca 



(16) 



where C is an arbitrary integration constant, which later 
on will be taken to be \/C — 1/{SttG), and Fx is the 
derivative of F{X) with respect to the kinetic term. All 
in all, there are three independent equations, namely the 
equation of the scalar field, the conservation equation for 
dust, and the second Friedmann equation. For a given 
Hubble parameter as a function of the red-shift H{z), 
the time evolution of the scalar field can be found using 
the set of the basic cosmological equations. This will be 
useful later on (see equations (29) and (30), as well as 
the definitions (27) and (28)). 

The comparison of a theoretical model against super- 
novae data relies on the minimization of 



V^ {Pthjzj) - Pobs{Zi)Y 



ot 



(17) 



with p the distance modulus, p = m — M, where M is the 
absolute magnitude and m is the apparent magnitude. 
The theoretical apparent magnitude pt^ as a function of 
the red-shift z = — 1 + ao/a is given by 



p{z) = 25 + blogio ( ^^) 



(18) 



where the luminosity distance d^ for a flat universe can 
be expressed as 



r 1 

dL{z)^c{l + z) dx 

Jo H(x) 



(19) 



and c is the speed of light, while H{z) is the Hubble 
parameter as a function of the redshift. 

We also exploit the CMB shift parameter R, since it is 
the least model dependent quantity extracted from the 
CMB power spectrum, i.e. it does not depend on the 
present value of the Hubble parameter Hq. The reduced 
distance R is written as 



Jo 



1089 



dz/H{z) 



(20) 



where J7„i is todays value of the matter normalized den- 
sity, and we use the CMB shift parameter value R = 
1.7 ± 0.03, as derived in |ll|. 

Independent geometrical probes are Baryon Acoustic 
Oscillations measurements. Acoustic oscillations in the 
photon-baryon plasma are imprinted in the matter dis- 
tribution. These BAOs have been detected in the spatial 
distribution of galaxies by the SDSS [12] at a redshift 
z ~ 0.35. The SDSS team reports its BAO measurement 
in terms of the A parameter 



A{z = 0.35) = Dv{z = 0.35) 



V ^m^O 



0.35c 
where the function Dv{z) is defined to be 

1/3 

Dviz) 
Da{z) 



dL{z)/{l + z) 



(21) 

(22) 
(23) 



and we have used the value ASDSS{z = 0.35) = 0.469 ± 
0.017. We shall assume that the scalar field Lagrangian, 
C = —F{X), is chosen so that the solution of the cosmo- 
logical equations for the Hubble parameter H{z) leads to 
the luminosity distance 



, , , c z(l + az) 

driz) = 

^ ' Ho l + (3z 



(24) 



Then the CMB shift parameter as well as the BAO A 
parameter, for a given matter density firm are certain 
functions of the two parameters a,/3, and if we require 
that R{a,(3) = 1.7 and A{a,(3) = 0.469 the two parame- 
ters can be determined. Taking Vim — 0.29 we find that 
the model is viable as long as the two free parameters 
take the values 



a = 1.008 
;3 = 0.318 



(25) 
(26) 



and furthermore the distance modulus p cannot be distin- 
guished from the corresponding quantity when the mat- 
ter is ignored, see Figure 1 below. 
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FIG. 1. Distance modulus versus red-shift for a) a model 
without matter, a = 1.246 and /? = 0.446 (dashed line), and 
b) our model including matter with fi,„ — 0.29, a — 1.008 
and P = 0.318 (solid hue). 
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FIG. 2. Time evolution /(r) in our model including matter 
with fl^ = 0.29, a = 1.008 and /3 = 0.318 (solid line), and 
the polynomial of second degree P2{t) = 0.2385r + 0.3297r^ 
that best fits the solution (dashed line). 



To find the time evolution of the scalar field, we intro- 
duce the dimensionless scalar field and time as follovifs 



/ 



Ho 

Hoit-to) 



(27) 
(28) 



where to is the age of the universe today, and we use the 
second Friedmann equation. The dimensionless time r 
as a function of the red-shift is given by 



r(z) 



dx 



{l + x)E{x) 



(29) 



while the dimensionless field / as a function of the red- 
shift is given by 



/(^)-/(0) = 



dx 



m. 



V2E'{x) 
V2{l+x)E{x) (l + a;)3 



(30) 
where E{z) = H{z)/Ho, and is given in terms of the 
assumed (dimensionless) luminosity distance Dl{z) — 
{Ho/c)dLiz) by 



^(^) - ( d-z 



Dl{z] 



Dl{z) = 



l + z 
z(l + az) 
l + Pz 



(31) 

(32) 



with a, P having the values determined before, and we 
have used that = \/2X. Finally, eliminating the red- 
shift we obtain / as a function of r, and the solution can 
be seen in Figure 2 below. In the same figure we also 
show the polynomial of second degree that best fits the 
solution /(r) in the same interval, which is the following 



P2(t) =0.2385r + 0.3297T2 



(33) 



We see that the solution exhibits a global minimum at 
T ~ —0.35, which corresponds to a red-shift z ~ 0.48. 
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FIG. 3. Deceleration parameter versus red-shift for a) the 
standard KCDM model with Q.m = 0.27 and f^A = 0.73 
(dashed line), and b) our model including matter with fi,„ = 
0.29, a = 1.008 and /3 = 0.318 (solid line). 



This is the value for which the two terms in (|30l) cancel 
eachother. 

Before ending our discussion, let us compare the decel- 
eration parameter q = —d/{aH^) of the standard KCDM 
model with that of our model. Using the definitions, the 
deceleration parameter as a function of the red-shift is 
given by 



q{z) = -l + (l + z) 



H'{z) 
H{z) 



(34) 



and can be seen in Figure 3, where for KCDM we have 
used the values Vim = 0.27, J^a = 0.73. In our model 
the universe enters into the acceleration era at z — 1.16, 
while in the standard cosmological model this takes place 
at later times. 

In summary, in the present work we have assumed 
that the current cosmic acceleration is driven by a scalar 
field, the Lagrangian of which is a function of the ki- 
netic term only, and that the luminosity distance is a 



given function of the red-shift. Upon comparison with 
Baryon Acoustic OsciUations (BAOs) and Cosmic Mi- 
crowave Background (CMB) data the parameters of the 
models are determined, and then the time evohition of 
the scalar field is obtained by integrating the cosmolog- 
ical equations numerically. We find that the solution is 
very different than the corresponding solution when the 
non-relativistic matter is ignored, and that the universe 
enters the acceleration era at larger red-shift compared 
to the standard ACDM model. 



APPENDIX 

Here we briefly discuss how a given function f{x) that 
vanishes at the origin can be fitted by a polynomial of 
second degree P2{x) = ax^ -I- 6a; in a given interval xi < 
X < X2- We require that the deviation of the function 
from the polynomial is as lower as it can be, or in other 
words the integral 



Iia,b) 



dx{.f{x)^P2{x)f 



(35) 



must be extremized. If we set the first derivatives of / 
with respect to a, h equal to zero 



dl{a,b) 

da 
dlja, h) 

db 



= 
= 



(36) 

(37) 



we obtain a two-by-two system 

ha + hb = 
I^a + hb - 

where we have defined 



Ml 

A*2 



(38) 
(39) 



Mn = 


/ dxx''f 

Jxi 

I-X2 


(x) 


(40) 


^7n ^ 


/ dxx"" 




(41) 



for n = 1,2 and m = 2,3,4. For the given interval and 
function the above integrals can be computed, and solv- 
ing the simple algebraic system we find the values of the 
coefficients a, b given in the text. 
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